Macromolecular theory of solvation and structure in mixtures of colloids and polymers 



M. Fuchs 1 * and K. S. Schweizer 2 
1 Department of Physics and Astronomy, The University of Edinburgh, JCMB King's Buildings, Mayfield Road, Edinburgh 

EH9 3JZ, United Kingdom 

2 Departments of Materials Science & Engineering and Chemistry and Materials Research Laboratory, University of Illinois, 

Urbana, Illinois 61801, USA 
(February 1, 2008) 

The structural and thermodynamic properties of mixtures of colloidal spheres and non-adsorbing 
polymer chains are studied within a novel general two-component macromolecular liquid state ap- 
proach applicable for all size asymmetry ratios. The dilute limits, when one of the components is at 
infinite dilution but the other concentrated, are presented and compared to field theory and models 
which replace polymer coils with spheres. Whereas the derived analytical results compare well, 
qualitatively and quantitatively, with mean-field scaling laws where available, important differences 
from "effective sphere" approaches are found for large polymer sizes or semi-dilute concentrations. 



I. INTRODUCTION 

Mixtures of dispersed spherical particles and non- 
adsorbing polymers may be viewed as a model system for 
a wide variety of materials encountered in food products, 
biological systems or technological applications. In these 
systems, the "depletion attraction" is always present be- 
cause it has a purely entropic and universal origin. Its 
consequences can most clearly be studied in mixtures of 
colloidal hard spheres and polymer chains made up of 
hard units where only the entropic consideration of the 
packing of particles restricted by steric or excluded vol- 
ume constraints enters. 

Because of the fundamental nature of the depletion 
attraction it has been studied theoretically since the 
pioneering work of Asakura and Oosawa Q and Vrij 
P). Moreover, its effect on the phase behaviour had 
been observed much earlier S. The phase diagram 
of colloid-polymer mixtures has been constructed using 
the Asakura-Oosawa pair potential in an effective one- 
component thermodynamic perturbation calculation [Q, 
within a two-component dilute polymer free volume ap- 
proach j5|, and also with computer simulations Q and 
0, the latter based on a specific model originating in 
Refs. ||. The Asakura-Oosawa model consists of re- 
placing the polymer coils with effective spheres which 
can freely interpenetrate each other but not the colloidal 
spheres. This model has been further treated by liquid 
state theory has been extended to non- homogeneous 
situations P, 1C], and to perturbatively include polymer 
non-ideality [11|. The forces it predicts for dilute rather 
large colloidal spheres have been measured directly fl^ , 
|l3| , and phase diagrams for colloidal spheres appreciably 
larger than the polymers have been obtained, and agree 
semi-quantitatively with theory . 

More detailed experiments on the colloidal correlations 



5|Jlj] |l9| ], measurements of the second virial coefficient 
20f , and quantitative tests of the phase diagrams for 



larger polymer sizes 

|j|miDlll!r|2§, havc ' however ' de_ 

tected polymer correlations which are not contained in 
the mentioned approaches. Also, for small spherical sur- 
factant micelles the dependence of phase separation 
on polymer size is opposite to predictions of colloid ap- 
proaches |Q, |f| . There are at least two reasons for these 
discrepancies. First, polymer coils can deform close to 
particles and can thus fit into void spaces more effectively 
than spheres. Second, for higher polymer concentrations 
the coils start to overlap and the relevant polymer cor- 
relation length crosses over from the coil radius to the 
size of a mesh in the formed transient network. Both ef- 
fects are important if the polymer coils are not negligibly 
small compared to the particles, and both are neglected 
in the described theoretical approaches. The effects have 
long been understood from field-theoretic approaches to 
polymers in the limit of dilute colloidal particles. The 
deformability of the polymer coils affects the depletion 
layer of polymer segments close to particles, the result- 
ing insertion free energies, and the induced colloid pair 
interactions p5|-p9fl . For semi-dilute polymer suspensions 
the depletion layer and the induced interactions were ob- 
tained for both large and small colloids (2?j [33). Yet, ex- 
cept for in a highly idealized mean field thermodynamic 
perturbation calculation to hard spheres by Schaink and 
Smit 0, polymer field-theoretical approaches have not 
been extended to finite colloid concentrations. 

Recently we proposed a macromolecular liquid state 
theory for mixtures of arbitrary polymer to colloid size 
ratios J39], which, although it is not rigorous for dilute 
systems, presents a viable and first principles approach 
for finite densities. It is unique in its applicability to all 
parameter ranges concerning densities and sizes, and is a 
macromolecular generalization pq-p8[ of the interaction 
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site description introduced by Chandler and Andersen 
for small rigid molecules [^9|. This generalization has 
proven rather successful for pure, especially dense poly- 
mer systems and polymer alloys. Some results for di- 
lute particle mixture systems have been obtained within 
a simplification of the approach |fl||4l|], and could ratio- 
nalize several surprising aspects of measured second virial 
coefficients of small proteins p0[ |. Also, light scattering 
measurements of the colloid liquid structure could be de- 
scribed semi-quantitatively over all length scales without 
adjustable parameters p5| ]. 

In the present manuscript we analyze in detail the low 
density limits of this macromolecular approach. The rea- 
sons are threefold. First, by looking at polymer solutions 
containing few colloidal particles it is possible to com- 
pare with exact field-theoretic results and thus to test 
the approach. Second, by considering dilute polymers in 
a hard sphere solvent it is possible to make contact with 
the previous Asakura-Oosawa type approaches. Third, 
in these limits fully analytical solutions of the non-linear 
integral equations description are possible and provide 
insight into the theory, which also applies to the higher 
concentration states. 

A conceptually new closure (approximation) for the 
direct correlation function describing the packing of 
polymers close to repulsive walls or around hard col- 
loidal spheres has been introduced in ]35| ] which en- 
tails a medium-ranged colloid-polymer segment effective 
interaction. To capture the key physics the latter is 
required within the polymer reference interaction site 
model (PRISM) approach since a preaveraging approx- 
imation for the single polymer chain form factor is em- 
ployed for tractability reasons. In inhomogeneous sys- 
tems, however, the single polymer form factor depends on 
the distance of the polymer chain from interfaces or in- 
homogeneities. Considering a fluid of random walk poly- 
mers, Gaussian intramolecular correlations apply. The 
number of intersections of a random walk with a plane 
scales as y/N, where N is the number of steps or polymer 
repeat units. Without rearrangements, the number of 
contacts with a repulsive wall would scale identically, as 
the excluded volume constraint could be satisfied by just 
mirror inverting the overlapping polymer strands. This 
result follows from PRISM with the most simple excluded 
volume closure (of Percus-Yevick form) jl(J. Close to 
the repulsive wall, however, the (Gaussian) intramolecu- 
lar polymer correlations differ from the ones in the bulk 
solution as translational entropy can be gained by re- 
ducing the number of contacts with the wall to 0(1), as 
required for recovery of the ideal gas equation of state 
from the wall virial theorem. In order to describe the in- 
homogeneous system with one homogeneous polymer in- 
tramolecular structure factor, the rearrangements close 
to a colloidal particle need to be captured by an effec- 
tive colloid-polymer interaction which extends across the 
range where the polymer segments rearrange. We pro- 



posed a molecular closure convoluting the local (bare) 
segmental steric repulsion with a Yukawa weight where 
the range, or non-locality length, called A, is determined 
from thermodynamic consistency considerations. On the 
segmental level the assumption of short ranged effective 
steric interactions entered into a Percus-Yevick style ap- 
proximation HEl. 

Thermodynamic consistency correlates the structure 
on local length scales with long wavelength fluctuations 
and is a familiar concept in liquid state theories [|42| . In 
the present case, the polymer chemical potential at infi- 
nite dilution, closely connected to the insertion free en- 
ergy for adding dilute polymers to a particle fluid, is used 
to implement consistency as it provides one of the sim- 
plest measures of the tendency of colloids and polymers 
to mix. Moreover, this quantity also determines the free 
volume which is one of the input quantities in the colloid 
theory most widely used for large colloid to polymer size 
ratios [|| 

The outline of this paper is as follows: In section II the 
model of colloid polymer-mixtures is presented. Section 
III describes the solution of the integral equations in the 
two low density limits of interest. The thermodynamic 
consistency equations are solved and discussed in section 
IV. Section V then presents the results and discussions 
for the structure and thermodynamics of dilute colloidal 
particles in a polymer solvent, while section VI describes 
the opposite case of dilute polymer chains immersed in 
a hard sphere fluid. Conclusions are in presented section 
VII, and three appendices contain technical material and 
a discussion of alternative closure approximations. 



II. MODEL 

The binary mixture shall be described by its (matrix 
of partial) structure factors Sij (q) , where the index i — 1 
indicates the polymer and 2 the colloid component. A 
small molecule solvent is treated as a background contin- 
uum and enters only implicitly via the interaction poten- 
tials for the polymers and colloids. In principle, all par- 
tial structure factors are experimentally measurable by 
(labelling and) scattering techniques. The total density 
fluctuations are decomposed into single molecule contri- 
butions, described by a (diagonal) intramolecular form 
factor, t2>ij(q) = LUi(q)6ij, and intermolecular correlations, 
hij(q), resulting in: 

S(q) = guj(q) + gh(q) g . (1) 

An obvious matrix notation is used. The diagonal ma- 
trix of densities, Qij = giSij, gives the number density 
of colloidal particles and polymer segments. The pair 
decomposable excluded volume or steric interaction pre- 
vents the particles/ segments from overlapping: 
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,(r< ~(<r i +<r j )) = 



(2) 



where 02 = cr c is the colloidal hard core diameter and 
<J\ = Cp is the excluded volume diameter of a sin- 
gle polymer repeat unit (segment). The intermolecu- 
lar pair correlation functions, gij(r), are trivially con- 
nected to the total intermolecular, hij, correlation func- 
tions, gij(r) — hij(r) + 1. Carets in Eq. ([[]) denote 
Fourier-transformed quantities. The total density fluc- 
tuations of the interacting fluid are decomposed into the 
single molecule fluctuations and in an interaction part via 
a generalized Ornstcin-Zcrnickc, or Chandler- Andersen, 
equation QH^: 



(3) 



In a pre-averaging approximation the single polymer- 
molecule density fluctuations 
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are taken to be known a priori. As the colloidal particle 
is assumed to be rigid and solid, it acts as a point scat- 
terer, uj c — 1, subject to the steric restriction, Eq. If 
the effective interaction potentials, the direct correlation 
functions, Cjj(r) in Eq. (^|), were taken to be the bare 
pair potentials, Cij(r) = —Vij(r)/(kBT), then Eqs. ([l]j|) 
would correspond to the RPA approximation which is one 
of the simplest liquid state approximations for (dense) 
polymeric and simple fluids. Typically one finds that 
RPA solutions violate the excluded volume condition, Eq. 
(||) . Integral equations approaches like PRISM go beyond 
the RPA as they enforce the no-overlap condition Eq. (||) 
rigorously and determine the direct correlation functions 
from self-consistency equations implementing the (phys- 
ically) motivated expectation that the Cij(r) are short 
ranged and vanish beyond a few particle diameters. For 
the colloidal hard sphere component this corresponds to 
the well established Percus-Yevick (PY) approximation 
il: 



3 (r > <r c ) = 



(5) 



This closure and the excluded volume constraint, Eq. 
(0), together with the site-site Ornstein-Zernicke equa- 
tion, Eq. (0), result in a coupling of density fluctuations 
at different wave vectors, thus leading to non-linear inte- 
gral equations with a much richer mathematical structure 
than the simple RPA approximation. For the model of 
hard spheres, the only thermodynamic parameter is the 
packing fraction, (f> c = ^g c a^. 

Detailed studies of the PRISM equations for ho- 
mopolymer solutions and melts (37],|38| have established 
that the polymer site-site direct correlation function to a 
good approximation decays to zero beyond the polymer 



repeat unit size so that a correspondingly simple closure 
can be enforced: 







(6) 



Thus, the interaction between polymer macromolecules 
is made up of pairwise site-site segmental interactions, 
which are given by a spherically symmetric effective po- 
tential, which follows from the excluded volume con- 
straint, Eq. (||). Attractive interactions beyond the 
"athermal" model studied here can be included [pj], 
Site averaged quantities are considered, and therefore 
specific chain-end effects are neglected. 

As discussed in the introduction, the effective colloid- 
polymer interaction extends beyond the range of imme- 
diate overlap. This arises because of, and allows to ac- 
commodate, the change of the polymer conformations 
close to colloidal particles. As the exact direct correla- 
tion function is not known we suggested |3fj] a simple 
one-parameter extension of the PY closure (m-PY): 



c s cp {q) 

1 + q 2 X 2 



with 



JL ) = 0, (7a) 



which enforces excluded volume on the local scale by fix- 
ing c s cp {r) from the excluded overlap condition, Eq. (||), 
and from the requirement of short-ranged segmental in- 
teractions. On physical grounds, because —UbT c s cp {r) 
describes interactions on the segmental scale, one expects 
c cp( r ) ^0 be negative (repulsive) and to exhibit rapid vari- 
ations (on the segmental length scale) and smoother ones 
connected with the colloid size. In real space the closure 
clearly implies a smearing of the segment-colloid interac- 
tions over the distance A: 



Cc P (r) = Jd*s -Lj j-1 



e-l r - s l/ A c s cp (s) 



(7b) 



where the PY closure for c s cp {r) can be viewed as de- 
scribing unconnected polymer segments, and the (nonlo- 
cal) conformational constraints on the segment packing 
( "chain connectivity" ) close to colloidal particles are cap- 
tured by the spatial convoulution. 

The m-PY closure contains an undetermined param- 
eter, the length A, which can be expected to vary non- 
trivially with the physical system parameters, like den- 
sities or size ratio. As it captures the rearrangements 
of the polymers strands close to a colloidal particle, its 
magnitude should be of the order or smaller than the 
polymer correlation length (i.e. radius of gyration, R g 
for dilute systems, blob diameter or mesh size for semi- 
dilute conditions) and/or the colloid size. Also, the poly- 
mer conformational changes, and hence A, will depend on 
the volume taken up by the colloidal spheres. In order 
to achieve a parameter free a priori description, ther- 
modynamic consistency shall be enforced te determine A 
uniquely. The implementation of this well known con- 
cept within liquid state theories starts from the observa- 
tion that the PY closure |Ba|, A = 0, leads to results for 
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the solution free energies obtained from the compress- 
ibility theorem, d ® dg F ex - — —k B Tcij(0), where F ox - is 
the excess free energy per unit volume, which compare 
favourably with field theoretic results where available 
fhi|| . Thus, the excess chemical potential for inserting 
polymers into a hard sphere fluid, where 5/ii = -^-F cx -, 
as obtained via the "compressibility" route provides a 
rather A-insensitive reference quantity since it empha- 
sizes long wavelenth correlations. Especially the limit for 
vanishing polymer concentration shall be discussed, 



\g P =0 



d{f c Nc ep (q = ) ef c )\ ttp =o, (8) 



where = l/(fceT), and the expression per molecule is 
given. An independent, more local route to the inser- 
tion free energy will lead to strongly A-dependent re- 
sults, thereby allowing a sensitive determination of A 
from equating both expressions. The approach of ther- 
modynamic integration introduced for RISM approaches 
by Chandler 0] shall be used as it connects the pair 
correlation functions on local distances to the thermody- 
namic properties. The variation of the free energy when 
turning on the interactions via the Mayer / function, 



AC) 

Jaf3 



A 1 ) 



is the 



: C e [0, 1] i-> shall be used, where f^ p 

physical one and f^i belongs to some known reference 
system; here the Greek indices run over the colloid and 
the polymer segment sites. As only excluded volume in- 
teractions are present, and because the limit of vanishing 
polymer segment size, <r p — > 0, is of interest, the ther- 
modynamic integration is chosen to take a mixture of 
polymers and colloidal point particles of free energy Fq 
to the true mixture by growing the colloidal particles, 



crP = (a c . From Ref. 



(3(F - F ) 



HQ c Qp<J c 



one then easily finds: 



(9) 



where the g^ are the ^-dependent (e.g. via the volume 
fraction of the colloid particles) pair correlation functions 
which are evaluated at the distances of closest approach. 
Equation (|]) expresses that the growing colloidal spheres 
have to push against the pressure of the surrounding sys- 
tem (polymers and colloids) which — in a virial theorem 
analogy — is given by the probability of contact on the 
surface. Immediately, one obtains a second independent 
result for the chemical potential of Eq. (||) which, as ar- 
gued, depends on A strongly as the packing of polymer 
segments close to the colloidal particles, g C p{\{<Jc + c p )), 
enters crucially: 



\g P =0 



Trg c a c N 



?9 { S 
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\g P =0 



\e P =o 



(10) 



Even though it would be desirable to obtain A for all 
concentration ranges, especially its form for low poly- 
mer concentrations is required. On the one hand, even 
small amounts of polymers added to colloidal systems can 
strongly affect the phase diagram and colloid structure. 
On the other hand, if A's dependence on the polymer 
parameters is known then scaling considerations allow 
reasonable extrapolations of A into the semidilute region 
as will be shown in section V. Therefore, the two expres- 
sions, Eqs. (PJlO|), for vanishing polymer concentrations 
will be used to obtain A in section IV. 



III. SOLUTION IN LOW DENSITY LIMITS 



The specified model of colloid-polymer mixtures cov- 
ers all polymer and colloid density regions, and distances 
from the polymer repeat unit size (ca. 5 A) up to the 
collective correlation length approaching phase separa- 
tion (some /zm). Since small amounts of non-adsorbing 
polymer can alter the structure and phase diagram of a 
colloidal system quite appreciably, simplification to con- 
sider rather low polymer concentrations is of initial inter- 
est. Under these conditions there opens up a mesoscopic 
window where polymer segments are well separated but 
polymer molecules overlap and interact strongly. Such 
solutions are called dilute or semi-dilute [Q, and treat- 
ment of this regime is most conveniently done by per- 
forming the "thread limit" |l5|,|4(|, where the size of a 
polymer segment and the corresponding statistical seg- 
ment size, Z p /\/l2, are taken to be negligibly small; 
Op oc l p — > 0. In order to retain polymer molecules 
with a finite radius of gyration R g , the number of re- 
peat units is increased beyond bounds, N — > oo, such 
that £q = IpN = Rg/2 = fixed. As has been shown 
by a rigorous solution of the PRISM integral equations 
in |4^| , intermolecular excluded volume remains active if 
in parallel the monomer density is increased, g p — > oo, 
such that the number of polymer molecules per coil vol- 
ume (~ Rg) stays finite: (p p — 27t^-£q — fixed. The 
reduced polymer concentration ip p differs only by a nu- 
merical factor from the often used polymer packing frac- 
tion 7] p = ^fjj-R 3 g = Qplep ~ (fp/0.53, where g*, is 
the density when polymer coils start to interpenetrate. 
The mathematical thread limit of the PRISM equations 
corresponds to a scaling law description of the dilutc- 
to-semidilute crossover of polymer solutions and can be 
compared to field theoretic scaling laws and results . 
In both cases only mesoscopic parameters, the polymer 
molecule density and coil size, enter and all microscopic 
parameters, like a p , l p and g p <Jp, drop out. The ef- 
fective polymer-polymer interaction becomes of the Ed- 



wards delta- function type, c pp (r) 



pp 



(0)(5(r), where the 
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intermolecular excluded volume parameter, c pp (0), fol- 
lows self-consistently from the no-overlap condition, Eq. 
(||). Besides its use for dilute and semi-dilute polymer 
solutions, experience also has shown when applying the 
thread limit outside its rigorous range of validity that 
it describes qualitatively adequately the spatially coarse- 
grained features of concentrated polymer solutions and 
melts || . 

Two limits of the scaling function of the single chain 
form factor are known in general and for Gaussian poly- 
mers simplify to: u(q — 0) = TV and uo(qR g 1) — > 
(qlp)~ 2 - Note that the self-scattering term, which is 
present in the full u>(q) of Eq. (||) for q<j p = 0(1), is not 
accessible in the thread limit. In order to keep simple and 
analytically tractable equations, the full intramolecular 
structure factor shall be approximated by the standard 
Pade interpolation between the two asymptotes: 



w(q) 



N 



(11) 



As the single-polymer structure factor is an input to our 
PRISM approach, the use of random walk statistics in 
Eq. ( |ll|) for repelling coils can be considered an addi- 
tional technical approximation in order to achieve ana- 
lytical results. In order to capture effects of the non- 
trivial intramolecular correlations ("swelling" and "self- 
avoiding- walk statistics" ) caused by intramolecular ex- 
cluded volume, the equations (|l]) to ( |To|) c ould be solved 
numerically with an appropriate u>(q) |4q| . Also, effects 
specific to semi-flexible polymers and arising from local 
chain rigidity are neglected in Eq. (|Tl]), but could be 
incorporated into numerical studies. 

In the following (and on the rhs of Eq. (pi])) dimension- 
less units shall be chosen by using the colloid diameter 
as unit of length, <r c = 1. Then, the length scale ratio, 
£o = Rg/(V2&c), the relative polymer concentration, tp p , 
and the colloid packing fraction, <fi c , are the only remain- 
ing physical parameters. Further notational simplifica- 



tion is provided by defining S cc = S cc /g c , S pp = -^r—S ppi 



4 P = c s cp<ri/lp, and w(q) = (l p /a c ) 2 uj(q). 

The limit of considering only (semi-) dilute polymer 
solutions does not eliminate the non-linearities of the in- 
tegral equations for the polymer and colloid structure. 
Insights into the physics described by the m-PY PRISM 
equations and their full solutions can be gained by reduc- 
ing one of the densities further to a dilute limit where at 
most pair-wise direct interactions of the diluted species 
can occur. This linearizes the equations in the correlation 
functions of the diluted species and thus, as the correla- 
tions of the majority component are known, simplifies 
the analysis. These limits will be studied in the follow- 
ing, where in section III, the m-PY closure parameter A 
still is kept arbitrary. 



A. Dilute Colloids 

In the limit <fi c — ► the equations simplify as the colloid 
particles do not alter the structure of the polymer fluid. 
The collective polymer structure factor for Gaussian in- 
tramolecular correlations within PRISM equals |B8) : 



e 



i + q 2 e 



(12) 



where the polymer correlation length crosses over from 
ss R g at high dilution to the blob-size or density screening 
length for concentrations within the semidilute regime: 



1 _ 1 1 



1 + 2ip p 



(13) 



Note, that the neglect of microscopic length scales cor- 
responds to the assumption that the polymer concentra- 
tions in Eq. ( [l3|) are far smaller than melt densities, 
where segments of different polymer chains start to pack 
densely (typically ~ 30 — 40% of melt density). 

The polymer segment profile close to a single colloidal 
particle and the resulting packing of two colloidal hard 
spheres is described by linear equations: 



h cp (q) = 



cl(g) 



Spp{q) 



h cc {q) = c cc (q) 



q 2 X 2~PP 

Vp 



and 



Cp 



(q) 



2<o l + q 2 X 2 



hcp(q) 



(14) 



Closed equations for the pair correlation functions follow 
from and Eqs. (||Q and (§|) - (@). Even though these 
equations can be solved by straightforward conversion to 
differential equations, a simplified Wiener-Hopf factor- 
ization of Eq. (|L4|) is used in Appendix A because it is 
close to the solution technique of the full equations for 
arbitrary <fi c and ip p to be published elsewhere, and can 
be presented in a more concise way. In Appendix A.l 
it is shown that the functions fij(r) = r(gij(r) — 1) sat- 
isfy simple differential equations. The one describing the 
polymer segment density profile close to a single sphere 



(l + \d r )(l + £d r )f cp (r) =0 



for r > 



(15) 



with initial conditions f' cp (j) = g' cp (j)/2 — 1 = — 1 and 

f'Jpih) = 

denotes a derivative. Also the direct correlation function 



i + A + 0/(^0; where a prime 



at zero wave vector is obtained from Eqs. (A3 A9): 



>(0) 



6N£ 



r (l + 6f +12f + 6A(l + 20' 



12A 2 (l + 2£)) 
(16) 



For dilute polymer solutions, (f„ — > 0, the quantity 
—Nc cp (0) is the cross second virial coefficient describ- 
ing the mutual excluded volume between a hard sphere 
and a Gaussian polymer coil. 
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For the pair correlation function describing the proba- 
bility of two isolated spheres in the polymer solution to 
be at a separation r, the differential equation is: 



(1 + Xd r )f cc (r) + ^(u + v\d r )f cp (r - 1) 



PY result of Ref. |4Q] is recovered. The virial coefficient, 
as it still contains the parameter A, will be discussed in 
section V, after the determination of A from thermody- 
namic consistency. 







for r > 1 



(17) 



with initial condition f cc (l) + 1 = g e c(l)- The unknown 
parameters u and v follow from Eqs. (A£,A12), and the 
initial value is given by the colloid pair contact value at 
infinite dilution, d> c — ► 0: 



flcc(l) = 1 + 
Vp (2g 2 -4A 3 (l+2£)+4Ag(l+2£)+A 2 (3+4g-4£ 2 )) 



(18) 



where corrections of the order C(e _1 / A ) were neglected. 
The contact probability for two colloids in principle can 
be calculated exactly within m-PY, but the approxima- 
tion to neglect corrections of 0{e~ 1 / x ) simplifies the ex- 
pressions greatly, and for realistic values, (see below) 
where A < (mostly <C) 0.31, introduces errors of at most 
a few percent. 

From Eq. (|l5| ) and the initial conditions follows that 
9cp{r) consists of a superposition of two Yukawa tails, 
A^ r ^/r and A A e _r/A /r, where A$ > and A\ < for 
A < £, and that it describes a monotonous increase from 
zero to unity for | < r < oo. Similarly, Eq. (0) shows 
that gcc(r) is a superposition of these two Yukawa tails 
(with constants and B\) and of an additional secular 
term C\e~ r / X . Nevertheless, numerically it is found to 
decrease monotonously for all parameters. 

Of interest is the second virial coefficient i?2 which fol- 
lows from the (colloid partial) compressibility and mea- 
sures the strength of the polymer induced pair potential 



B. Dilute Polymers 



In the limit of adding only a few polymers to a dense 
fluid of hard spheres, the m-PY PRISM equations sim- 
plify to: 



h cp (q) 



1 + q 2 X 2 
h P p(q) = u' 2 q c pp {q) 



S cc (q) 



and 



2 - , n , He uc (q) t 

—hc P {q) , 



7T 1 + q 2 X 2 



(21) 



where the structure factor of the pure hard sphere fluid, 
S cc , is given in PY-approximation. The unperturbed 
length scale characterizing the polymer correlations is 
given by £o as for an isolated polymer, but colloid me- 
diated interactions arise as the polymer coil is forced to 
squeeze into the voids between the spheres. In Appendix 
A. 2 it is shown that the functions fij(r) = r(gij(r) — 1) 
satisfy linear integro-differential equations with known 
integration kernels. The average polymer segment den- 
sity around colloidal spheres increases according to: 

,1/2 

(1 + Xd r )(l + Zod r ) f cp (r) = I2<j> c $ / ds q cp (s)f cc (r - s) 

J -1/2 
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where the initial conditions are: f' cp {\) = 
-1 and 
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where the result for hard spheres, B^ s = expresses 
a purely steric repulsion, and negative values of Bi indi- 
cate the presence of a net attractive effective interaction 
due to "depletion" . One finds from the solution of Eq. 

(0) 



2A£o(l-</>c) 2 



Here, the hard sphere functi on, f cc (f) is known or can be 
easily determined from Eq. ( A21 ) in Appendix A. 2. The 
new factor function is given by: 



q cp (r 
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2f (1+0-4A 5 (l + 20 + A 4 (l + 20 2 
+ 4A£ 2 (l + 3£ + 2£ 2 ) +4A 3 (l + 4£ + 6£ 2 
+ 2A 2 £ (3 + ll£+12f 2 + 4£ 3 )] , 

where again corrections of order C(e _1//A ) were neglected 
but could be evaluated straightforwardly. For A = 0, the 



and vanishes outside this range; the parameters are given 
in Eq. (A25). Again, the direct correlation function at 
zero wave vector, which is needed for the thermodynamic 
^ ^calculation, can be obtained in explicit form (see Ap- 
^ pendix A.2): 



(20) 
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The microscopic colloid-segment hard core interaction 
causes a repulsive spike at contact which leads to an os- 
cillatory behaviour for very large wave vectors: 



km -> °°) 



-2ttCOs(-) 



C + 2(A + £ O )(1 + 20 C ) 



(1-0 C ) 2 



(26) 



This remnant of the local segmental interactions also is 
present at finite polymer concentrations, and for C — > 
is given by Eq. with the replacement £o — ► 

£. The polymer pair correlation function follows from 
(|A27|jA29|jA3l|) for r > 



(1 + \d r )(l + £d r )f pp (r) = 2 & e- r /f° + z x e- r / x 

+ 120 c £ o 2 fK ds q cp (s)(r - s)(g cp (\r - a\) - 1) , (27) 

where the initial conditions are f PP {0) — and f pp (0) = 
— 1, and the parameters are given in Eq. jA30| ). Clearly, 
by means of the last term the colloidal spheres imprint 
their local packing structure onto the dilute polymers. 
In order to explicitly calculate the pair correlation func- 
tions, numerical integrations of Eqs. (^,^) are most 
convenient using the known hard-sphere function, which 
can also be obtained from Eq. (A21). 

Also of interest is the intcrmolccular excluded volume 
parameter which describes a colloid modified effective in- 
teraction for polymer segments on different chains: 

(o)£o _ b £o+ z >.+ z to _ l 

-SttZJ A+£ (1-0) 4 (A+Co) 



(1 



'■to (1-0(1-36))) 



+ 3A 2 0(l + 2£ o -0 (1-4&)) +A (1-0) 
x (1-0 (2-9£o -12£ 2 ) +0 2 (l-9£ + 24£ 2 ))] . 

It simplifies to c pp — Sirl*/^ for the pure polymer sys- 
tem, and is connected to the polymer second virial coef- 
ficient via: 



B P 2 = -~ h pp (q = 0) = -\N* ( c pp (Q) + ^ (c cp (0) r S 



hard sphere fluids. The second, Eq. jic|), uses very local 
information captured in the contact values of the pair cor- 
relation functions. Within integral equation theories for 
the structure of many-body systems both results gener- 
ally do not coincide, and in the m-PY closure this aspect 
of "thermodynamic consistency" is used to determine the 
effective interactions length scale A from equating both 
results. 

The ideal gas result for adding point polymers, £o — > 0, 
immediately follows from Eqs. (^,^5|) and the observation 
that A oc £o (corresponding to the polymer appearing as 
an inpenetrable small sphere to the colliod) has to hold 

in this limit: Nf35fi p c) | ep=0 = - In (1 - 4> c ) + O(£ ). This 
change in translational entropy when dissolving a point 
polymer in a sphere solvent of packing fraction C is con- 
nected to the free volume fraction accessible to the poly- 
mer, e~ Nl35 ^ = Vf /V = 1 — C . The free volume for 
finite colloid but vanishing polymer concentrations is a 
central quantity in the phase studies of the free volume 
approach by Lekkerkerker et al. |p[. In order to arrive at 
the same ideal gas limit from the polymer-colloid contact 
probability in Eq. @, since ^g cc (l)| ep=0 = #(£o) as 
follows from Eq. (|l8|), the following result has to hold 
for the contact probability of polymer segments with the 
colloidal surface: 



.1 
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for Co -> 



(30) 



(29) 



which becomes B\ 
are present. 



~N 2 c pp (0) oc Rg when no colloids 



The contact value should be microscopically small, van- 
ishing in the mesoscopic limit, N — -> oo. In the thread 
m-PY PRISM limit described in section II, where micro- 
scopic parameters were scaled away, the colloid-polymer 
pair correlation function agrees with this result. Thread 
PRISM correctly predicts on mesoscopic scales that 
(28^ cp (I) — > 0. In the case of (semi-) dilute polymer so- 
lutions a scaling connection of the segment-segment con- 
tact value to the mesoscopic scaling law description could 
be shown, which enabled one to estimate the microscopi- 
cally small contact value from an evaluation of the meso- 
scopic pair correlation function at a microscopic separa- 
tion pS| . From the parabolic polymer segment density 
profile close to a colloidal particle and Eq. j23|), one 
realizes that within m-PY this again is possible. The 
mesoscopic pair correlation function (i.e. the thread g cp ) 
evaluated at a microscopic separation, a p , because of 
9cp{\ + Op/2) = 5<?cp(i) *5p> predicts the correct scal- 
ing behaviour of the contact value, Eq. (63), 



IV. ENFORCEMENT OF THERMODYNAMIC 
CONSISTENCY 

The solutions for the structure of colloid-polymer mix- 
tures open up two very different calculation routes to the 
thermodynamic properties. One, via the long wavelength 
fluctuations, termed "compressibility route" , is suggested 
in Eq. (ph for the insertion free energy of polymers into 



9c P ( 1+ o ap ) 



1 ol 



1 



2 A£o 1 - 



for Co -» , (31) 



and thus recovers the ideal gas result for the chemical po- 
tential from the "wall virial" route. Evidently, however, 
the microscopic distance could only be calculated if the 
full PRISM equations also including microscopic length 
scales were solved. As this has not been achieved yet, in 
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Ref. 1351 for simplicity the distance a 



cr? = 21^ was chosen, 
because then the distance over which point like polymers 
rearrange close to a colloid sphere becomes numerically 
identical to the polymer correlation length: A — £o for 
£o ^ 0. Note that this holds for all colloid packing frac- 
tions, and that a p has to be chosen only once. 

It is the major difference of the novel m-PY closure 
to the previously studied PY description, where A = 0, 
and (as can be expected) to prior numerical solutions of 
the PRISM-PY equations @|§, that thermodynamic 
consistency in Sfi p is possible in a scaling sense. For 
A = in the PY-case, the mesoscopic pair correlation 
function predicts a scaling of the colloid-polymer con- 



tact value as g cp { 



^) = 5p/(£o(l - <j> c )) for £ -> 0, 
thus violating Eq. ( p0[ ) and the ideal gas insertion 
free energy. Obviously, the PY closure overestimates 
the segment density of polymers, even of an isolated 
polymer, close to walls or colloidal spheres, as a fac- 
tor >/~N too high contact probability is predicted. This 
comparison justifies the motivation and the physical in- 
terpretation of the m-PY closure in section II. Addi- 
tionally, the PY-approximation A = predicts a seg- 
ment density profile which increases much more strongly, 



^H(r-i)+0((r-i) 2 ),com- 
i"\2 1 mil™ i\3^ 



(l-0c)^ 



pared to <? c °(r) t (r - \) 2 + 0((r - from 

Eq. (|22|), and thus strongly underestimates the polymer 
induced depletion attraction. This follows since for small 
separation of two colloidal particles, PY predicts a linear 
increase of the polymer density in the gap, whereas the 
correct ^5[Q and m-PY result for Gaussian polymers 
is a quadratic increase. 

Extension of the expansion for small £0 leads to the 
higher order result A = £o~ 4£q+£qA(0 c ) + . . ., where the 
dependence on the sphere fluid packing density indicates 
an anomalous behaviour of A for 4> c — > 1; i.e. one finds 
A(0) = 16 but A(4> c ) -► - (1 ^ )2 for C — * 1. As A is re- 
stricted to positive values, a different expansion scheme is 
required for high colloid packing fractions, and for large 
polymers or small colloidal particles. It is described in 
Appendix B. The calculation leads to A — > A(0 C ) for 
£0 — * 00 where: 



a- 1 = r x 



1 + 2</> e Ai 

1 — 6r Or 



(33) 



where Ai = (v5 + 1), £ — £0 for vanishing polymer con- 
centration, and for convenience dimensional units are re- 
stored. This Pade approximation satisfies the thermo- 
dynamic consistency condition from Eqs. (| 10) up to 



relative errors of 15 % for all parameter values, see Fig. 
[jj. For small polymers, A follows the polymer correlation 
length. For large single polymer coils, A is expected to 
be determined by the pure hard sphere fluid correlation 
length beyond which density fluctuations are screened, 
which is c -dependent. The result in Eq. (33) can be 
shown to be in excellent agreement with this intuitive 
idea. 

The expression in Eq. ( |33| ) immediately suggests an 
extrapolation for A from its present calculation at van- 
ishing polymer concentration to finite tp p . If the poly- 
mer correlation length of the dilute situation, £ , is re_ 
placed with the full density dependent one, £ from Eq. 
(|l3|), then A also is known for 4> c = but now arbi- 
trary polymer concentrations, tp p > 0. This replacement 
is suggested as the role of the single polymer molecule 
structure factor, to(q), in Eq. ( pi] ) is taken over by the 
collective one, S pp (q), if the densities are changed accord- 
ingly. This standard procedure from polymer-scaling ap- 
proaches thus substitutes the blob-size or density screen- 
ing length in place of the chain size for semi-dilute situ- 
ations |Q In section VI, it will be shown that this ex- 
trapolation achieves thermodynamic consistency for the 
insertion free energy of adding colloids to a polymer so- 
lution, i.e. for quite a different thermodynamic quan- 
tity as originally considered for the determination of A. 
In Ref. |p5| , this way of extrapolating A to finite poly- 
mer and colloid concentrations was suggested and used 
with £ = £(</> c , f pi £,o) being the full polymer correlation 
length, which depends on <p c , tp p and £o- ln the two cases 
considered here it simplifies to £(</> c , = 0, £0) = £0 and 
£(0 C = 0, ^,£0) = £o/(l + 2<P P ), see Eq. @. 



A(0 C 



0.80. 



1/3 



T (1-0 C ) + O.78...(1-0 C ) 2 + 



(32) 



which verifies the interpretation of A as the range of poly- 
mer segment rearrangement, because it is always less 
than the smaller of the two molecular sizes, and becomes 
vanishingly small if no free volume is available for the 
polymers. Exact determination of A from Eqs. (p|]To|) for 
all parameter values is difficult, because thermodynamic 
integrations are required. Thus we proposed the approx- 
imation to interpolate between the known exact limits 



V. RESULTS FOR COLLOIDS DILUTED IN A 
' POLYMER SOLVENT 

If only a small amount of colloidal hard spheres is dis- 
solved in a polymer solution, then the structure of the 
polymer fluid is not affected. The intramolecular density 
fluctuations or the form factors for the individual poly- 
mers, (jj(g), are rather well understood from field the- 
oretic considerations |5^]. PRISM describes the pack- 
ing, i.e. the intermolccular correlations, of the polymeric 
macromolecules which are simplified to Gaussian chain 
molecules in the present work. The intermolecular pair 
correlation function exhibits the well-known "correlation 
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hole", which shows that polymer molecules for entropic 
reasons softly repel each other. The correlation hole has 
a non-trivial structure on the length scale of the size of 
the molecule (radius of gyration) and also for shorter dis- 
tances, on the mesh- or blob scale characterised by the 
density screening length £. 

Close to the diluted colloidal particle the polymer seg- 
ment density is less than in bulk resulting in a depletion 
layer of varying width. Its form is described by Eq. Jl5| ) 
and shown in Fig. |^. The width is given by the polymer 
correlation length £ as long as the particle is much larger. 
Only for dilute polymer solutions, this correlation length 
agrees with the polymer size, £o- I n the semidilute re- 
gion, it is given by the blob size or mesh width. If the 
particle becomes smaller than the polymer correlation 
length, then the depletion layer width crosses over to the 
particle diameter. This is shown in the inset in Fig. |^. 
If very small particles are immersed in the polymer solu- 
tion, then the depletion layer has an additional power law 
tail due to chain connectivity correlations which becomes 
g cp (r > 1) = 1 - (| + 2? for £ - oo ||. 

In the dilute polymer solution limit, the PRISM re- 
sult compares rather well with field theoretic calculations 
from Ref . |^5) . While the width is of the same order in 
both results, especially the small distance power-law in- 
crease of the density profile, g cp (r) = \g" p {\){ r — ^f) 2 + 
. . ., can be compared as it obeys a scaling law, where the 
(universal) amplitude of PRISM (semi-) quantitatively 
agrees with the known field theoretic (FT) limits: 



-k (PRISM) 



2 A 



h (PRISM) 



£ (FT) £ 
£ (FT) £o - 







oo 



(34) 



The PRISM result for the polymer profile in Eq. (|lf 
extends into the semi-dilute concentration region where 
mean field calculations for the case of R g <§; <J C are avail- 
able p0| |. Again, the amplitude, g" p (^) — > for £ — > 0, 
agrees verifying its universality (2^, . Only in the case 
of very large polymers, m-PY overestimates the deple- 
tion effect by ca. 20% in Eq. (S3). Figure [2] shows that 
the width of the depletion layer for semidilute solutions 
is set by the polymer correlation length (for £ <C er c ), 
and that the molecular size, R g plays no role for the de- 
pletion layer; this also is evident from Eq. (^) where 
only the correlation length of the collective polymer fluc- 
tuations appear. Using an effective pair potential ap- 
proach ^j, |B| determined under dilute polymer condi- 
tions strongly overestimates the range of the depletion 
interactions in semidilute polymer solutions. 

The polymer induced depletion attraction becomes ap- 
parent from the colloid pair correlation function describ- 
ing the probability of two isolated colloidal spheres to be 
at a distance r, Eqs. (JTJ,|l^). As seen in Fig. ||, for close 



distances and especially at contact, this probability is in- 
creased above the random value of unity. The g C c(f) de- 
crease from the contact values, Eq. (|l8|), monotonously 
without oscillatory features or layering. Thus, the ef- 
fective or induced potential, (3V cfi -(r) = — \ng cc (r), is 
attractive and monotonic, and does not exhibit repulsive 
barriers. The contact value scales like g cc (l) — 1 ~ ^ 
in either limit of very small or very large polymers, i.e. 
£o — > (oo). The distance characterising the decay of 
g C c{f) is closely connected to the width of the depletion 
layer as seen in the inset of Fig. |^, and both depend on 
the polymer parameters only via the nonlocality length 
A. The correlation between two colloidal particles which 
are much smaller than the polymer mesh width falls off 

for large distances as g cc (r > 1) -► l + (l + 2/Ai) 2 7r ' p ^ gp 
for £ — * oo, indicating a weak but long-ranged attraction 
p7[ . It becomes screened at the polymer or blob size 
only. 

A measure of the tendency of colloidal spheres to dis- 
solve in polymeric fluids is given by the chemical poten- 
tials or solvation energies. From the compressibility route 
in analogy to Eq. (ih, one finds: 
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(c)| 



dg' p c cp (q = 0, £>p)| ec=0 ■ (35) 



From the free energy (see Eq. (Q)) required to grow the 
colloid particle from a point to its actual size one finds: 



( 6 ) U=Q = ttc^ f^ CC 2 ff (C)»A U=0 . (36) 



P6fi. 



For dilute polymer solutions, where £ = £oj the result 
from both routes obeys thermodynamic consistency (up 
to errors of 10 %) and agrees with field theoretic results 
up to 20 % as shown in fig. ||. From the local packing 
information, Eq. (pq), an explicit result can be found: 
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(37) 



which has a number of polymer specific features. For a 
fluid of small polymers or far into the semi-dilute regime, 
£ — * 0, inserting a colloidal sphere costs the free energy 
of creating its volume by doing work against the osmotic 
pressure of the polymer fluid. For larger polymers, the 
form of the polymeric coil enters, and in the limit of very 
large polymer chains the chemical potential becomes in- 
dependent of R g and scales linearly with the length over 
which polymer segments need to be rearranged, which is 
just the colloid diameter a c . As the added sphere sees 
local strands of the polymer network only, the result for 
small spherical colloids becomes independent of the poly- 
mer's size (R g ) or degree of polymerization, and also in- 
dependent of the mesh or blob size. Thus it becomes 
independent of the polymer molecule concentration p9| : 



P5/i c \g c= o — > 7r0 p 0vJpAi for £ — > oo 



(38) 
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which holds in the dilute and semi-dilute region. This 
result also quantitatively compares favourably to the 
known behaviour in the dilute limit 0S(x c \ gc =o — > 
4irQ p a c lp. For the semidilute region, the two routes in 

Eqs. (HH) predict /3<5^ c) | ec=0 | /35/i s) | ec=0 and 



e c =o 



pn. = -c p ipl 



for 



(39) 



which explicitly states the connection to the osmotic 
pressure, II, of the pure polymer system, which is given 
in PRISM approximation for Gaussian polymers in Eq. 
( |39| ) . Note, that the extrapolation of A to finite polymer 
concentration leads to thermodynamic consistency also 
in these cases (as claimed at the end of section IV), as 
both routes, Eqs. (|35|,|36), lead to Eq. J3g|) identically, 
and the results in Eq. (39|) only differ by a numerical 
prefactor, which could be absorbed into a density depen- 
dent microscopic length, a p ((p p ). 

The mean polymer induced depletion attraction for 
two colloidal spheres can be quantified using the second 
virial coefficient which is accessible experimentally. 
Representative results are shown in Fig [| The m-PY re- 
sult, Eq. (p9|) , expresses that the effective colloidal pair 
potential depends on the polymer concentration and the 
polymer to colloid size ratio. For small polymers, Eq. 
( |l9| ) simplifies to: 



B c 2 



1 - 



12<p p 3 (28 Ai - 53) ip p & 



so 



l + 2(p p 4(l + 2<p P Y 



for £ 
(40) 



which indicates an appreciable attraction between the 
two colloids. Note that the third term is positive and 
predicts a weakening of the attraction when the poly- 
mer size starts to grow. Increasing the polymer density 
strengthens the attraction which saturates at a finite neg- 
ative value. The reliability of these results is discussed 
at the end of this section and in appendix C where alter- 
native closures are examined. 

For very large polymers, again an effective two-colloid 
attraction is observed. 



B c 2 



1 



3(2 + Ai) 2 ip p ^ /a c ) 
2\\ (1 + 2ip p ) 2 

(24+18A!+3A 2 V P 
2X1 (l + 2p„) 



for £ 



(41) 



As seen in Fig. |^, two hard spheres immersed in a 
fluid of much larger polymers feel an induced attraction 
which varies non-monotonically with polymer concentra- 
tion and increases with polymer size as previously found 
based on the A = PY closure p0[ . Its origin is purely 
entropic as the added particles hinder conformational 
fluctuations of a polymer strand they are embedded in. 



For large distances, the conformational entropy loss at 
both particles is additive. On shorter distances, how- 
ever, chain connectivity restricts the conformational fluc- 
tuations of a polymer molecule even without added parti- 
cles. Thus the additional entropy loss is smaller if the two 
particles are close. This induces an effective long-ranged 
attraction among the colloidal particles. For small con- 
centrations, adding polymer strengthens the effective at- 
traction. Around the dilute-semidilute crossover density, 
however, the effective range of the polymer induced at- 
traction crosses over from the polymer size to the mesh 
size, and thus starts to decrease appreciably. The latter 
effect dominates in the semidilute concentration region 
and therefore B% decreases again. It is important to men- 
tion that a mean field approach [^3| employing a RPA 
approximation, where £ ~ would be too crude to 

handle this competition and would miss the increase of 
£?2 above the overlap concentration. Apparently, a sim- 
ple superposition approximation which decouples the de- 
pletion layers around each of the colloid spheres misses 
the long-ranged, but weak, attraction |Q. The mini- 
mum value of B2 for the considered Gaussian polymer 

statistics is — (£o/o"c) ~ — 0.17(i? s /i? c ), which is 

asymptotically deeper than the value —0.50(R g / R C ) 0A0 
for self-avoiding walk polymer statistics |2i|. The more 
swollen polymer molecules in good solvents apparently 
are more open to the particles and allow stronger in- 
terpenetration so that the induced colloid attraction is 
smaller 

7 Deep in the semidilute density region, cp p — > 00, the 
second virial coefficient again saturates, B^/B^ 5, — ► —5, 
indicating a finite effect of the pair- wise attraction which 
is independent of the size ratio. This limit corresponds 
to a vanishingly small polymer blob size. Yet, for large 
polymers there can open up a window, where adding 
polymer only slightly changes B^, in a manner which 
is much weaker than at the overlap concentration but 
still differs from the asymptotic value. This is the poly- 
mer concentration range where the polymer mesh width 
is not yet negligible compared to the particle size, and 
there a broad maximum develops: 



Bl 
Bf s 



f° r fp ~* 00 1 Co — > 00 with £ = const. 

(42) 



The scaling function b% is shown in the inset of Fig. 
^. For semidilute polymer concentrations, the deple- 
tion attraction among two colloidal spheres depends non- 
monotonously on the ratio of the blob to the sphere size. 
There is an optimal blob size which roughly equals the 
sphere radius, where the induced attraction is minimal, 
or B% maximal. This makes physical sense since for 
a c /2 « £ the particles can "just fit in" the polymer mesh 
spaces without distorting it. For larger polymer corre- 
lation lengths, the crossover to Eq. (ET1) sets in because 
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the range of attraction increases. Unexpectedly, however, 
also for smaller £, becomes more negative, presumably 
because the depth of the attraction increases. 

There does exist a caveat for the above results on the 
colloid-colloid interactions 

f[ §,11,111. Liquid state the- 
ory (and also [Q) with PY-closure, density functional ap- 
proaches to the colloid structure, and free energy based 
approaches Q , all underestimate the depletion attraction 
in situations where it far exceeds fcsT. This error, which 
presumably affects the results for dilute colloidal parti- 
cles in a dense solvent of small polymers, </> c — > 0, £o <^ °c 
and if p 3> 1, arises from an inherent linearization of the 
depletion potential in the considered binary-mixture ap- 
proach. Therefore it is fortunate that a thermodynamic 
consistency condition can be formulated explicitly ad- 
dressing the accuracy with which the depletion attraction 
is handled. There exists another independent expression 
for E>2 which follows from Eqs. (p|Jl^) and the definition 
Eq. ©: 



Hg(/o lrf cc 2 ^S ) "(c/2) 



(43) 



f?c=0 



This result follows from the local packing information 
and shows a very different, and unphysical, scaling com- 
pared with the result from large wave length fluctuations 
for small polymers and/or high polymer concentrations: 

^ 2 Hg — > ■^('^p/S.o) 3 ■ This throws severe doubts on the 
result in Eq. (^) and on our treatment of the deple- 
tion attraction in these two cases, £o — ► or tp p — ► oo, 
where from a pair potential point of view the depletion at- 
traction should be arbitrarily large. Thus m-PY PRISM 
with the PY closure for colloid-colloid correlations ap- 
parently cannot capture the strong induced attraction 
in these two limits of theoretical interest as has been 
pointed out previously |3^,^0|. As the second virial co- 
efficient can be considered a worst case example of this 
failure, appendix C examines this issue using different 
theoretical approaches. Reassuringly, the results for the 
colloidal second virial coefficient for large polymers are 
recovered semi-quantitatively from the two routes within 
m-PY PRISM. The scaling in the limit of very large poly- 
mers is recovered exactly, i.e. the term diverging linearly 
with £ in Eq. (f4l|). Equation ([l3|) also predicts that 
a scaling law exists for the second virial coefficient of 
colloidal particles in a polymer mesh which semiquanti- 
tatively compares with the one of Eq. (p^), see Fig. [|. 
Intriguingly it also shows the non-montonic behaviour of 
B%(0 f° r intermediate £. 



VI. RESULTS FOR POLYMERS DILUTED IN A 
HARD SPHERE FLUID 



While numerous field theoretic results for polymer so- 
lutions exist, including for the question of dissolving one 
or two particles in dilute solutions, and could serve for 
tests of the m-PY results in the previous section, to our 
knowledge little is known about the packing of dilute 
polymer chains into dense particle fluids. The case of 
present interest is when the particles are much larger than 
the repeat units of the polymer chain, and both com- 
ponents are immersed in a small molecule background 
solvent treated as a continuum. If the amount of added 
non-absorbing polymer is small, then the structure of the 
colloidal fluid is not changed. In the present treatment 
it is given by the reliable and easily improved PY-thcory 
description. 

Whereas the segmental depletion layer of a single poly- 
mer (or a semidilute solution) around a single (spherical) 
particle exhibits a monotonic dependence on the distance 
to the surface of the particle, adding further particles 
forces the polymer to squeeze into the open spaces. Thus 
the probability of finding polymer segments at a distance 
r from the center of a colloidal particle, g C p{r), develops 
an oscillatory structure whose period is correlated with 
the colloidal size. In Fig. |] the evolving layering of the 
segment density is shown. For high particle densities, the 
polymers pack tightly into the voids and thus are close to 
the particles. The depletion layer, even though present in 
principle, is restricted by the colloid spacing, and varies 
with external parameters. Thus, the assumption of an 
effective pair potential becomes inappropriate, because 
it requires that the range of the polymer induced colloid- 
colloid potential, which naturally is connected to the de- 
pletion layer, itself becomes particle density dependent. 
The inset of Fig. [6] shows that there is a strong correlation 
of the colloid-polymer interaction range, A, to the width 
in the depletion layer. While varying the polymer size by 
4 orders of magnitude, neither the width nor A change as 
strongly. Moreover, for somewhat higher colloid solvent 
concentrations the width becomes a unique function of 
A, which only mildly splays out if <p c is decreased. Both 
quantities furthermore arrest at finite values in the limit 
of £o — > oo, since the relevant length scale is then the 
colloid size. 

If the solvation free energy for adding non-adsorbing 
polymers to a fluid of spheres is considered, Eqs. fljiol), 
then the ideal gas limit result discussed in section IV 
only holds for point polymers, £o — > 0. As expected and 
shown in Fig. ^, it becomes more difficult to add larger 
polymers because less free volume is available for them. 
The compressibility expression, Eq. (g) can be integrated 
analytically yielding: 
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(44) 



which was compared to the virial route in Fig. || For 
larger polymer coils, the added macromolecule looses 
conformational entropy when squeezing into the fluid in- 
terstitials and thus the free energy cost increases. Yet, 
its increase is not as rapid as if the polymer was a sphere 
because the polymer chain can rearrange. This becomes 
especially important for large polymers because they can 
wrap around the fluid particles. Therefore, if Eq. (Q) is 
compared with the corresponding result from scaled par- 
ticle or PY-theory for the addition of spheres with size 
R g to a sphere fluid || then two qualitative differences 
appear for large R g 3> o c . First, only a quadratic scaling 
with the polymer size results, 
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whereas for (large) added spheres the chemical potential 
scales with the volume of the added species. The be- 
haviour in Eq. (filf ) is connected to the large-£ scaling 
of the chemical potential for adding spheres to a poly- 
mer solution, discussed in Eqs. ( |3^j39| ), because both 
quantities are determined from the polymers' ability to 
deform around a particle. In Eq. ( |45| ) this explains why 
the rhs becomes independent of the degree of polymeriza- 
tion and linearly dependent on the colloid size. Second, 
the increase with particle fluid density is weaker than 
the corresponding sphere mixture result, which would 
predict a cubic divergence for cf> c — ► 1. Thus the dif- 
ference of both approaches becomes more important at 
higher sphere concentrations. Note that the free volume 
expression used Ref. || is connected to Sfj, p \ gp —o via 
a = Vf/V = e r Nl3& ^^=° . The input to this approach, 
the chemical potential for adding a sphere of radius R g to 
a hard sphere solution, is included in Fig. |7[ and for large 
polymers strongly overestimates the free energy cost for 
insertion. 

The pressure exerted by the particle fluid on the poly- 
mer also manifests itself in the dilute limit intermolec- 
ular packing of the segments of two polymers which is 
described by g pp (r) for ip p — > 0; examples are shown in 
Fig. @. The potential-of-mean force between segments 
is given by — fc^Tln g pp (r), while its polymer molecule 
analog is of order N 2 times larger. The slight repulsion 
of polymers in solution which causes the correlation hole 
is overcome at the distance connected to the depletion 
layer in g cp ; as seen by comparison with Fig ^|. For larger 
distances the polymer segments are pushed together and 
therefore pack more densely than random. This corre- 
sponds to a strong tendency of polymers to cluster and in- 
terpenetrate even though g p <C Q p *- For small polymers 
an oscillatory pattern develops connected to the distance 



of voids between the particles. For larger polymers, this 
segmental layering is smeared out and after an initial rise 
g PP {r) monotonically decays to its random value. In the 
limit where the particles are much smaller than the two 
polymeric molecules, but of course still much larger than 
the segmental repeat unit size, the two polymers entwine 
strongly. For £ ^> 1 and ip p — > 0: 



b P (r) - 1 + (/(0c) - 1) 

where f(6 ) - ^(6A 1+ i-40 e ) 
wnere J{q> c ) - 2 (i-0 c )(i+20 c ) > 



(46) 



follows from Eqs. (£2|,^7|). Within the coil radius, for 
r < Co, the distribution of segments from two chains, 
is almost constant and does not exhibit the self-similar 
power-law behaviour, 1/r, of open (Gaussian) fractals. 
Only for distances larger than the coil size, r > £o, does 
the intermolecular segment distribution decay to uncor- 
rected packing. The probability of segments of different 
polymer chains to be close, as measured by the pair cor- 
relation function contact value foo(4> c )i is determined by 
the particle packing fraction, and for <f> c > 0.11 becomes 
(much) higher than the bulk density. This increase is 
a precursor of the demixing transition at finite polymer 
concentrations pq] . For <p c w 0.5, the density of segments 
from other chains within a R g distance of a tagged seg- 
ment is more than an order of magnitude greater than in 
the absence of colloids. This may have significant con- 
sequences for intermolecular processes such a chemical 
reactions or energy transfer between polymers added at 
dilute levels to colloidal suspensions or porous materials. 

The intermolecular excluded volume parameter, de- 
— Cpp(0), or the related polymer 



fined as Vp Xcl -(<p c ) — — <- PP 
molecule second virial coefficient 



i?2, determine the im- 



portance of the excluded volume interaction on the poly- 
mer packing on local and macromolecular distances, re- 
spectively. Results are shown in Fig. [| On local dis- 
tances, it appears plausible that a polymer chain should 
experience the identical steric repulsion from a segment of 
its own backbone as from another of the chemically iden- 
tical macromolecules. Thus, the result in Eq. ( ^8| ) can be 
used to discuss the effective excluded volume parameter 
induced by the depletant particles. For vanishing colloid 
concentration, where Vp Xc1 ' and are directly related, 
the interpretation of polymer chains as repulsive spheres 
is recovered in PRISM because Vp Xcl -((j) c ) oc l/R g and 
therefore B\ oc Rg. Note, that our present model thus 
does not describe O-solvents, where B\ = 0, even though 
a Gaussian intramolecular structure was assumed a pri- 
ori, but rather describes polymers in athermal solutions 
with the technical simplification to treat the polymer 
chains as random walks. For results with self- avoiding 
statistics, see, e. g. Refs. p§, [36]. The result in Eq. 
(Eq), which to a very good approximation simplifies to 
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indicates that the effective local segmental interaction 
stays repulsive for all size ratios and particle fluid con- 
centrations. Asymptotically for small depletant particles, 
it becomes independent of the macromolecular size, and 
saturates to a positive value which increases with (f> c . The 
depletion attraction induced by the polymer on the par- 
ticles overcomes the effectiveness of the particles to push 
polymer segments together on a local scale. This is found 
even though, mesoscopically, two polymers are induced to 
intertwine, and for somewhat higher polymer concentra- 
tion, phase separation into a polymer gas and polymer 
fluid phase sets in |35|. On macromolecular distances, 
however, the effective polymer-polymer molecule inter- 
action can become very attractive if the particle density 
is high, or the colloidal particles are large: 
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For small polymers, the result known from PY theory for 
hard sphere mixtures is recovvered. It is negative and 
describes the tendency of the polymer point particles to 
cluster. Increasing the polymer size somewhat increases 
this tendency. For large polymers, however, a finite col- 
loid density is required in order for the colloid induced 
attraction to overwhelm the segmental repulsion. The 
virial cocficient becomes negative only for <fi > 0.108 . . . 
for £ > 1- 



VII. CONCLUSIONS 

We have studied athermal colloid and polymer mix- 
tures and considered the dilute limit of one species using 
a binary mixture approach which treats the hard spheres 
and Gaussian polymer coils on an equal footing. The 
macromolecular liquid state theory uniquely addresses 
the structural correlations over the wide range of length 
scales from the polymer repeat unit size to the molecu- 
lar sizes. Packing of hard spheres is handled using the 
reliable Percus-Yevick approximation. Polymer specific 
effects are captured which arise from the ability of the 
polymer coils to deform close to and around the parti- 
cles. Also, the appropriate polymer correlation length 
appears in the description because the formation of a 
polymer mesh is captured. The polymer specific effects 
become important as soon as the size of the coils is not 
negligible relative to the colloid radius. 



The comparison of the m-PY approximation with rig- 
orous field theoretic results in the dilute colloid limit 
presented in section V serves to validate our approach 
which offers the unique possibility of extension to higher 
densities. Interesting and novel predictions arise for the 
colloid induced pair interaction of dilute polymers which 
are much larger than the particles. Strong interpenetra- 
tion is predicted (see Fig. ||) finally leading to fluid-fluid 
phase separation fl35| . 

In this context it is important to stress that the ap- 
proximation of Gaussian single chain correlations is done 
for purely technical reasons in order to achieve analytical 
results. Of course this entails that the majority of the 
scaling predictions derived in the present and previous 
work j35) bears the wrong exponents if applied to poly- 
mers in a good solvent. Intramolecular excluded volume, 
thermal attractions and a self-consistent determination 
of oj(q) will be included in future numerical studies, and 
in a number of cases the corrected scaling predictions 
have already been pointed out ||^, Q . 

Colloidal dispersions containing added free polymer 
are often described by integrating out the polymer de- 
grees of freedom in order to derive an effective colloidal 
Hamiltonian Q]. However, for interacting polymers there 
exists no small parameter, and the induced many-particle 
interactions among the colloids do not generally termi- 
nate at a pair-wise description, and the polymer den- 
sity dependence of the parameters in the effective inter- 
actions can not be neglected. Both assumptions enter 
the description of polymer-colloid mixtures using effec- 
tive pair-potentials like the Asakura-Osawa model 
and are justified only for low concentrations of polymers 
much smaller than the colloid particles, and preferrably 
dissolved in O-solvents which mimimize the excluded vol- 
ume interaction. The discussion of semi-dilute polymer 
solutions in section V clearly identifies the importance of 
the density dependent polymer blob size or mesh-width, 
which also appears in O-solvents (44). The importance 
of induced many-body interactions, which naturally are 
contained in the presented binary mixture approach, was 
discussed in Refs. |6[^5|. Treatments of binary mix- 
tures of spheres have been used to describe the colloid- 
polymer mixtures, replacing the polymer coils by spheres 
of equal or similar size and neglecting their direct inter- 
actions ("phantom sphere models") p,H]. These capture 
induced many-body interactions but neglect the deforma- 
bility of polymers around particles and the crossover of 
the relevant polymer length scale from the coil size to 
the blob size when reaching semi-dilute concentrations. 
Whereas these effects had been well appreciated when 
adding dilute particles to polymer solutions, their study 
in concentrated particle solutions in the present approach 
clearly reveals their importance whenever the size of the 
polymer coils is not negligible. Their description in re- 
cent "soft-colloid" approches to polymers is yet unclear 
p7|,|5q]. An interesting quantity in this context is the 



13 



free energy cost for adding polymers to hard sphere flu- 
ids, because this is the central input for the widely used 
free volume theory of Lekkerkerker et al. || . As Fig. 
shows polymer deformability leads to strong deviations 
compared to results based on replacing polymer coils by 
hard spheres. 

A central ingredient of the present approach consists 
in enforcing thermodynamic consistency for the polymer 
insertion chemical potential S/i p \ gp= o. The length scale 
A over which polymer segments are allowed to rearrange 
close to the particles is determined from equating expres- 
sions for 8 [ip from long wavelength and from local packing 
information. Although the latter compares favourably 
with exact results for dilute systems, quantitatively we 
consider the former coarse-grained approach ( "compress- 
ibility route" ) to be more reliable, and also preferentially 
use it for other thermodynamic quantities within the 
present approach. Our reasoning rests on three observa- 
tions: (i) The used scaling law approach (thread model) 
to semi-dilute polymer solutions is coarse-grained a pri- 
ori so a zero wave vector thermodynamic route is natu- 
ral. Recently it has been connected to a self-consistent 
Gaussian field theory ]59[| . (ii) The compressibility route 
leads to thermodynamic results far less dependent on the 
specific closure as can be seen by comparing the new m- 
PY with previous PY results 0|(j|,|4l]]. (Hi) In a mi- 
croscopic calculation the universality, in field-theoretic 
sense, of the long-wavelength structure could be shown 
explicitly pTfl , in contrast to the local structure route 
where microscopic parameters remained as prefactors. In 
the present context this implies that a full PRISM calcu- 
lation could lead to a local matching length in Eq. j3l| ) 
dependent on the excluded volume size. Thermodynamic 
consistency, nevertheless, should prove a useful concept 
also for the description of other systems, like polyelec- 
trolyte mixtures [p2|| , where PRISM-based numerical ap- 
proaches have already given interesting results. 
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APPENDIX A: FACTORIZATION OF THE 
INTEGRAL EQUATIONS 

Three-dimensional Fourier transformations leading to 
f(q) of functions f(r) depending on the radius r only, 
shall be simplified to one-dimensional ones by: 



f(q) = 2%l dre' qr f(r) where f(r) = dssf(s) 

J\r\ 

(Al) 



1. The limit </> c -> 

Inserting Eq. (|l2|), the left equation of Eq. ( [ilj ) be- 
comes: 



k2 = e 
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(A2) 



(1 + AV)(1 + £V) ' 
which suggests an ansatz for c s of the form: 

c 12 (q) = 2ire- lq / 2 (u a + i\qv a ) + (1 + iq\)(l + iq£,)qi 2 (q) 

(A3) 

This shifts the problem of finding cf 2 (r) for — i < r < \ , 
to the problem of finding the (Wiener-Hopf factor) func- 
tion (712(f), where 



gi2(<?)=27r/ dre lqr q 12 (r) 



(A4) 



and qi 2 (r) = elsewhere. From the required symmetry, 



(q) = c 12 (— q), and Eq. (A3), follows the continuity of 



qi2 at the upper boundary: 



.1 



912(2) = 



(A5) 



Inserting Eq. (|AJ) into Eq. (|A2j) leads to 
(l-i«A)(l- iq£)h 12 (q) -e q 12 (q) 



2 . 2Kee- lq/2 {u a + i\qv a ) 



(l + iq\)(l+iqt) 
(A6) 

where by simple inspection the right hand side has no 
pole in the lower complex g-plane, ?sq < 0. Thus for 
r > — i the Fourier-back transform of Eq. ( [A 6] ) , remem- 
bering Eq. (Al), can be performed and leads to 



(1 + Xd r )(l + id r )h 12 {r) = e qi2 (r) for r > 



1 

2 ' 

(A7) 



For r > |, this leads to Eq. ( |l5| ) and in the overlap 
region, the excluded volume condition gives: 



r + A + £ = ey ia (r) 



1 1 

tor < r < — , 

2 2 



(A8) 



which can easily be integrated using Eq. (A5) and pro- 
vides the initial condition for the derivatives of f\ 2 at 
r = i mentioned in the main text. Also the zero wave 
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vector value of the direct correlation function in Eq. Jl6| ) 
follows from Eq. (A3) and the symmetry requirement on 
c , which fix: 



922 (0) 
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(A16) 



e +2A|(l+e) + X 2 (1 + 2Q 

e 



(A9) 



The ansatz Eq. ( A15 ) is useful because it cancels a 
number of the poles in the lower complex g-half plane in 
Eq. (All). The required short range of the direct cor- 
relation function, c 22 (r > 1) = 0, however, at first sight 
is violated by Eq. ( A15| ), which appears to indicate an 



This also predicts a Dirac-delta spike repulsion situated 
at the surface of the excluded region, which results in the 
large-g behavior: 



exponential tail, rc 22 (r) = const. e~' r '/ A for \r\ > j. Re- 
quiring the constant to be zero, i.e. requiring the residues 
of the rhs of Eq. (A15) at q — ±i/X to vanish, fixes the 
initial value of the factor function: 
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Inserting the ansatz Eq. ( |A3| ) at — q into the second 
equation of Eq. (|l4| ) leads to: 

t . h 12 (q)<p p ,2Tre lq/2 (u a - iXqv a ) (1 - iq£,)qi 2 (- 

"22 — C 2 2 — ^— {- 
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2< 



1 + AV 



(1 + iqX) 
(AH) 



where Eq. (A6) shows that the last term on the rhs has 
no pole for < except for a single pole at q = —j. 
Introducing the constants 



X(u a + Vg) 

2(A + 



and v = — |- z 
A 



(A12) 



also the first term on the rhs of Eq. ( All ) can be decom- 
posed into poles and zeros in different half-planes because 
of: 



(u a - iXqv a ) (1 - i£q) u - iXqv 

— z — 

1 + iXq 1 — iqX 



1 + AV 
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On the lhs of Eq. ( All ), as c 22 (r) vanishes for r > 1, its 
factorization can be achieved using a constant w and an 
undetermined function q 22 (r) which vanishes outside the 
range < r < 1, and whose Fourier transform is given 
by: 



q 22 (q) = 2tt dr e^ r q 22 {r) 
Jo 



(A14) 



Respecting that c 22 is a symmetric function of q, the fol- 
lowing ansatz is required: 



- / \ _ 2 TTW 1 <?22(q) , <?22(-g) 

c 22 {q) - 1+A 2 g 2 + jzixj + 1+lXq 

<P P !j 2 (■}i2(?)+27rze-"'/ 2 )(gi2(-g)+27rze"'/ 2 ) 
27Tfo 1 + AV 
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where the requirement of absence of a term of the form 
1/q 2 for q — > 00, which would correspond to a divergence 
C22(f — ► 0) ~ 1/r, determines the constant w from the 
value of q 22 at r = 0: 



+ 2^e 1 /( 2A )gi2(y ) + 2 7 rze- 1 /( 2A >g 12 (^)) , (A17) 

which is a linear equation in 922(0). Unfortunately, the 
||i)Ll expression for 922(0) turns out rather complicated 
and only simplifies if the fact that A <C 1 holds is recalled 
in order to neglect corrections of the order 0{e~ x / x ). 

The factorization of C22 expressed in Eq. (A15) simpli- 
fies the Fourier-back transformation of Eq. (All), as the 
poles in the lower complex g-plane were identified in Eqs. 



(A13,A15). Thus for r > one finds after Fourier-back 
transforming that 

(1 + Xd r )h 22 (r) + ^(u + vXd r )h 12 (r - 1) = q 22 (r) , 

(A18) 

holds, which proves Eq. ( |l7|) and, together with Eq. 
(A17), also proves Eq. (|18|), because of 322(1) = 
?22(1)/A. In the overlap region, it leads to: 

<l 22 ( r ) = r + X + ^-(ur - - + vX) for < r < 1 . 



ft 



(A19) 



The changes in the contact value upon adding more col- 
loidal particles, i.e. (dg cp /dg c )\g c =o, can be determined 
in an explicit form from the Wiener- Hopf factorization of 
the full non-linear equations |6(J . 



2. The limit ip p -> 

For (f p — > 0, the colloid structure factor agrees with 
the PY-solution for hard spheres, which can be written 
in terms of Baxter's factorization function as: 

- -1 

S 22 = (I - q 22 (-q)){l - q 22 (q)) where 

q 22 (q) = 2tt f dr q 22 (r) e^ r 
Jo 

r 1 a 

= 2tt dr{-{r 2 -l) + B{r-l))e lqr , (A20) 
Jo 2 
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-3<f> c /2 
£2(i_0 c )2. 

Correspondingly the pair correlation function satisfies: 



with the coefficients: A — g 2 (i-^fe and B 



Considering the explicit expression for /ii 2 (<?) from Eqs. 
(P, |A20| , |A22| ): 



/ 22 (r) -2w ds q 2 2(s)f22(r - s) = 
Jo 



for r > 1 . 
(A21) 



The first of Eq. (21) can be simplified with the ansatz: 



ci 2 (<?) = (i + ig£o)(i + i( i^)Qi2(q) 

+ e-^/ 2 (l - <7 22 ( (Z ))2^( Ub + igAu 6 ) , (A22) 

as from Baxter's solution it is known that (1 — g 22 (g))(l + 
£ 2 /i 22 (q)) has no poles for < 0. The function gi 2 (r) 
again is assumed to vanish outside the overlap region 
such that Eq. (A.4) holds. The required symmetry 
of ci 2 (q) and the known properties of g 22 further show 



that Eq. 



) holds and that vb = £oQi2(— h) an d 



Ub = -2irq 2 2(0)Xvb + v b {l + 2A + A/£ ). Inserting Eq. 
( A22j ) into Eq. (|l]) and closing the Fourier-integrals in 



the lower complex- q half plane thus leads for r > 

(l + £oW + A$.)Mr) 
= $<H2(r) + 12^ 2 f\ ds qi2{s)h 22 {r - s) . (A23) 



This leads to Eq. ( gg ) for r > | as / y (r) = —d r hij(r) 
for r > 0. Within the overlap region the equation deter- 
mining qi2 follows: 

4 1 

r + & + A = tfq[ 2 (r) + 12^ 2 /" ds q 12 (s)(r - s) , 



(A24) 
with pa- 



which is solved by the result given in Eq. (£ 
rameters: 

l-<f> e (1-6 A- 6£ ) 
a = ~ — 5 
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From the discontinuity of q' 12 at r = |, the initial condi- 
tion Eq. (p3|) results. Solving for the parameters 
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one can obtain the Eqs. ( p5| , p6| ) from Eq. (A22). The 
second of Eq. ( pl| ) can be rewritten upon insertion of Eq. 
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one recognises that the first term of Eq. ( A28) multiplied 
with (1 — (722(9)) contributes no poles in the lower g-half 
plane. Thus the poles for < 0, which are the only 
ones contributing for r > 0, from the second part of Eq. 
(A27) can be identified: 



120 c ^ e - i "/ 2 (l-g 22 (g))( Mi ,+ i :<;A^) 
(l+iqX)[l+iqio) 
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I ^ where the parameters are given by 
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( p7| ) follows by transformation and differentiation. The 
last unknown quantity, cn(q = 0), can be obtained from 
Eq. (A27) because the rhs vanishes as cos 2 (q/2)/q 8 , as 
can be seen from Eqs. (|21 A23). Therefore, the discon- 
tinuity of the third derivative at r = of the transform 
of WqCn needs to be balanced by Cn: 



9 r 3 /i 1 i(r)| r= o = -9 2 / 11 (r)| r=0 
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which together with Eq. (g7|) fixes the polymer excluded 
volume parameter as given in Eq. (pq) 

The change in the colloid-contact value upon adding 
polymers, {dg cc /dg p )\g p= o 1 can be obtained from the 
equations at finite density and leads to a rather unwieldy 
expression ]6C| ]. For this work, only its limits are re- 
quired: (dg C c/dg p )\g p= o = 0(£o) f° r small £ , and ne- 
glecting terms 0(A 2 ) 
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APPENDIX B: STEPS IN THE 
DETERMINATION OF A 

In the limit £o 1 one can use the physically mo- 
tivated expectation that A becomes small compared to 
£o and is of the order of the colloid diameter. With the 
ansatz, A = A(c/> c ) for £o — » oo, the expressions in Eqs. 
(E3LE5|) simplify. With Eq. ( A32 ) , which entails a further 



expansion in A < 1 up to linear order, the contact route 
expression Eq. ( ^0| ) becomes: 



l+2x 



An ordinary differential equation for A(0 C ) is obtained 
by considering: 



~ afloat; ~ |)^p c) le P =o] = , 



ratios £o or higher polymer concentrations) because the 
linearization in Eq. (CI) underestimates the strong at- 
traction in these cases. For larger polymers, however, in 
agreement with the consideration of thermodynamic con- 
sistency in section V, the PY linearization of Eq. (CI) 
is qualitatively and even quantitatively appropriate. The 
minimum in B^ at the overlap concentration, as well as 
the semi-dilute scaling law, are present in the HNC (and 
A = PY |40) results also. The closer quantitative agree- 
ment of the HNC results for B^ with the compressibility- 
route results within m-PY is one of the prime reasons 
we favour this route Eq. (p0|), over the virial type one, 
Eq. (fl3|). Moreover, whereas the latter was found to be 
(Bl\mphysical in the PY closure for the colloid-polymer in- 
teraction, A = the former which rests upon 
c C p{c[ = 0), show, identical qualitative behaviour as the 
m-PY. This can easily be observed by specifying the gen- 
eral result Eq. ( filf ) to A = 0. The novel m-PY closure 
predicts quantitatively stronger depletion effects than the 
(B2) an d moreover correctly describes the depletion layer. 



whose asymptotic solutions are given in Eq. (B2I) 



APPENDIX C: COMPARISON OF FROM 
OTHER CLOSURES 

The failure of thermodynamic consistency in the col- 
loid second virial coefficient B^ for situations where the 
effective colloid pair potential is strong, warrants a com- 
parison with results obtained with other closures. Be- 
cause of previous polymer-colloid mixture work in Refs. 
JlO|,|d| we consider two alternative approximations: the 
PY closure for the colloid-polymer interaction, which cor- 
responds to A = 0, and the HNC closure for the col- 
loid correlations. Study of the latter is motivated by the 
known deficiencies of the PY closure for mixtures of hard 
spheres, which are (partially) corrected by the HNC clo- 
sure 

The result for the colloid pair correlation function us- 
ing the HNC closure can easily be obtained from the m- 
PY result via Eol: 



7 HNC/. 



OO-i 



for r > <r c 



(CI) 



Because the quantity (3W(q) — QpC^. p (q)S pp {q) in Eq. 
( |l4| ) can be identified as (Fourier transform of) the poly- 
mer induced potential [p2| , the HNC closure agrees with 
the exact virial expression for colloidal spheres interact- 
ing with this (effective) pair potential and the a c /R g — > 
oo blob scaling mean field results |30|. This holds, be- 
cause j3W(r) = — .9cc _PY ( r ) f° r r > a c as follows from 
Eq. (14). Here, ,9cc~ PY ( r ) i s the pair correlation of the 
m-PY approach presented in the main text. From Fig. 
|lO| one deduces that, as discussed in section V, the m- 
PY results are not reliable for a c ^> £ (i.e. small size 
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FIG. 1. Polymer molecule excess chemical potentials, 5fi p \ ep =o, as functions of the size ratio £o from the two different routes 
used to enforce thermodynamic consistency, for three colloid packing fractions as labelled. The solid lines give the results from 
long wavelengths, Eq. (B), the dashed lines the corresponding ones from local packing, Eq. (|lo|). The inset shows the relative 
errors, which for all parameters are smaller than 15 %. 




FIG. 2. Polymer-colloid pair correlation function, g cp (r), exhibiting the polymer segment depletion layer close to an isolated 
colloidal sphere for various polymer correlation lengths £, as labelled. The inset shows a double logarithmic plot of the width 
w of the depletion layer, defined by g cp {r — | + w) = i, as a function of £; thin dashed lines mark the asymptotes w — > 1.66£ 
(w — > 1.1) for £ — > (£ — > oo) respectively. 
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FIG. 3. Colloid pair correlation function, g C c(r), for two isolated hard spheres in a solution of polymers with size ratio £o = 5 
for the denoted polymer concentration cp p . The inset shows as a solid line the range w of the depletion attraction estimated 
by g C c{r — 1 + w) = i(l + g cc (l)) as a function of the colloid-polymer interaction length A. Curves for various polymer 
concentrations (ip p = 0.1, 1, 2 ,10) overlap, while the size ratio £o runs between 0.01 < £o < 100. The dashed line compares 
the width of the depletion layer in g cp (r) from Fig. 0. 




FIG. 4. Colloid excess chemical potentials at infinite dilution, pS/j, c [ fle =a, as functions of the polymer to colloid size ratio £o 
for four polymer densities as labelled. The solid lines give the results from long wavelengths, Eq. (|35[), the dashed line is the 
result from local packing, Eqs. (^,^), in the dilute polymer limit. The short dashed line compares the field theoretic result 
known in the dilute limit [E5[ . The inset shows the relative errors in thermodynamic consistency for the four polymer densities. 
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FIG. 5. Colloid second virial coefficient B% in units of the hard sphere result, B^ s = as a function of th e poly mer 
concentration for various size ratios £o as labelled. The solid lines present the results from long wavelengths, Eqs. ( |2o| , ^o|j4ll ) , 
the dashed lines present the results from local packing, Eq. (^3|). The curves are cut in order to prevent overcrowding the figure. 
The inset shows the semidilute scaling law, 62(C) from Eq. (^), which applies to situations where the polymer correlation 
length £ (blob or mesh size) is of the order of the colloid diameter. The solid line again follows from Eq. (EOh and the dashed 
one from Eq. (fisj). 




FIG. 6. Polymer-colloid pair correlation function, g cp (r), for a polymer added to a hard sphere solution; colloid packing 
fraction <j> c and size ratio £0 as labelled. The inset shows the width w of the depletion layer (defined as g cp {r = \ + w) = ■=) 
versus the interaction range A for the labelled colloid packing fractions cj> c . The curves are parametrized by the size ratio £0 
which varies as 0.01 < £0 < 100 with increasing w and A. 
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FIG. 7. Polymer excess chemical potentials at infinite dilution, <5/x p | ep =o (in units of l/(N/3)) as functions of the colloid 
packing fraction for four size ratios £o as labelled. The solid lines give the results from long wavelengths, Eq. (^), the dashed 
lines the corresponding ones from local packing, Eq. The short dashed line compares the equivalent PY result for inserting 

a sphere of radius equal to R g as used in the phantom sphere approach of Ref . H . 
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FIG. 8. Polymer-polymer pair correlation function, g pp (r), for dilute polymer molecules immersed in a hard sphere solution; 
the colloid packing fractions <f> c and size ratios £o are the same as in Fig. Hand as labelled (<f> c — 0.5 for £o = 0.1 and 1). The 
thin lines present the asymptote for large polymers, Eq. (^6|), evaluated for £o = 5 and 1 at <fi c = 0.5, while the inset shows the 
intermolecular polymer segment contact value f(4> c ), which determines this asymptote. 
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FIG. 9. Normalized polymer second virial coefficient _&f* = -Bf/(47r£o) as a function of the colloid packing fraction 
for various size ratios £o as labelled. The inset shows the normalized polymer intermolecular excluded volume parameter, 
«p = «p Xcl '£o/(87rZp), as a function of the polymer size, £o, and the (indistinguishable) approximation Eq. (^7]); for the inset, 
the colloid packing fraction is <f> c = 0.2. 




FIG. 10. Colloid second virial coefficient B% in units of as a function of the polymer concentration for various size ratios 
£o as labelled. The solid lines are the compressibility results from Fig. ^ whereas the dashed lines give the corresponding HNC 
results. The inset presents the latter in a larger window, showing their drop to large negative values for polymer concentrations 
far above the overlap one. 
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